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ABSTRACT 

We present MERLIN maps of the 22-GHz H 2 O masers around four low-mass late- 
type stars (IK Tau U Ori, RT Vir and U Her), made with an angular resolution of ~15 
milliarcsec and a velocity resolution of 0.1 km s -1 . The H 2 O masers are found in thick 
expanding shells with inner radii ^6 to 16 au and outer radii four times larger. The 
expansion velocity increases radially through the H 2 O maser regions, with logarithmic 
velocity gradients of 0.5-0.9. IK Tau and RT Vir have well-filled H 2 O maser shells with 
a spatial offset between the near and far sides of the shell, which suggests that the 
masers are distributed in oblate spheroids inclined to the line of sight. U Ori and U 
Her have elongated poorly-filled shells with indications that the masers at the inner 
edge have been compressed by shocks; these stars also show OH maser flares. MERLIN 
resolves individual maser clouds, which have diameters of 2 - 4 au and filling factors 
of only ^0.01 with respect to the whole H 2 O maser shells. The CSE velocity structure 
gives additional evidence the maser clouds are density bounded. Masing clouds can be 
identified over a similar timescale to their sound crossing time (~ 2 yr) but not longer. 
The sizes and observed lifetimes of these clouds are an order of magnitude smaller than 
those around red supergiants, similar to the ratio of low-mass:high-mass stellar masses 
and sizes. This suggests that cloud size is determined by stellar properties, not local 
physical phenomena in the wind. 

Key words: masers - stars: AGB - circumstellar matter - stars: kinematics - stars: 
mass-loss - stars: evolution. 


1 INTRODUCTION 

Stars which start life with a mass M* of one or a few Mg 
usually end up as 0.6 - 1.0 Mg white dwarves (WDs). Most 
of their mass loss occurs during and at the end of the Asymp¬ 
totic Giant Branch (AGB) stage, at rates of M ~ (10~ 8 
10 -4 ) Mg yr -1 , resulting in dusty, thick circumstellar en¬ 
velopes (CSEs). Mass loss on the AGB determines the fi¬ 
nal stages of evolution of low- and intermediate-mass stars 
and significantly contributes to the chemical evolution of 
galaxies, in particular providing up to 80 percent of the dust 
(Whittet 1992). 

AGB stars contain a degenerate WD-like core sur- 
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rounded by a shell enriched with CNO-cycle elements and 
other products of high-temperature nucleosynthesis, which 
are brought to the stellar surface by convection and major 
dredge-up events known as thermal pulses. As stars enter 
the AGB their optical pulsation periods P become longer 
and more regular, of the order of a year. The outer layers 
are cool (typically 2000 - 2500 K) and very extended (ra¬ 
dius R t <1 1 au), with a luminosity L* of a few 10 3 Lg. 
The CSE becomes thicker during the AGB lifetime of a few 
times 10 4 — 10 5 yr (Jura 1993). More massive stars have 
shorter lifetimes but lose mass at a higher rate (Vassiliadis 
& Wood 1993, 1994; Blocker 1995). Finally, the remnants of 
the stellar shell heat up and are lost in a superwind, leading 
to a planetary nebula (PN) surrounding a WD. 

van der Veen & Habing (1988) located Miras, OH/IR 
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Table 1. The properties of the stars and their circumstellar envelopes. 


Source 

(i) 

Position (J2000) 

(hh mm ss.sss+dd mm ss.ss) 
(2) 

V, 

(km s -1 ) 
(3) 

P 

(days) 

(4) 

my 

(magnitude) 

(5) 

d 

(pc) 

(6) 

Mco 

M 0 yr" 1 

(7) 

Pm ax 

(km s -1 ) 
(8) 

r CO 

(au) 

(9) 

IK Tau 

03 53 28.83 +11 24 20.5 

+34.0 

470 

10.8 - 16.5 

265±| 5 

2.6xl0“ 6 

21.3 

6320, 26000 

U Ori 

05 55 49.169 +20 10 30.69 

-39.5 

368 

4.8 - 13.0 

260 ± 50 

2.3xl0 -7 

7.5 

7000 

RT Vir 

13 02 37.982 +05 11 08.38 

+18.2 

155, 

7.4 - 8.7 

133 ± 20 

1.3xl0“ 7 

8.8 

(2035 x 1100) 




112 ± 11 





±320 p.a. 30° 




& 170 ± 7 






U Her 

16 25 47.471 +18 53 32.87 

-14.5 

406 

6.4 - 13.4 

190 ± 70, 

3.4xl0 -7 

11.5 

5400 






385 ± 5 
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NOTES: 

(2) The position of IK Tau is accurate to better than 1", the positions of the other stars are accurate to ~ 10 mas. 

(3) V'lsk measurements derived from OH and SiO maser data, accurate to better than 0.5 km s —1 . 

(4) RT Vir varies irregularly. OH data (E01) suggests multiple periodicity. The periods of the other stars are accurate to within 10 days. 

(5) IK Tau and RT Vir are Zodiacal and have poorly sampled optical periods; IK Tau is also heavily obscured by its thick CSE. 

(6) In this analysis we adopt a distance of 133 pc for RT Vir and 266 pc for the other sources. 

(7) We used the distances given in note (6) and the most recent CO data available to calculate Mqo by the method of Groenewegen 
et al. (1999). Almost all available CO and IR data from the literature give results which agree to within a factor of 5. 

(8) Maximum CO or thermal SiO expansion velocity references as for (6). The discrepancies between these values and previous 
measurements is ss Tg km s —1 . 

(9) Maximum radius of CO emission, adjusted to our adopted distances. 
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(2001); G98 Kholopov et al. (1998); H97 Hipparcos Catalogue, The (1997); K85 Knapp & Morris (1985) K87 Kirrane (1987); K98 
Knapp et al. (1998); K99 Kerschbaum & Olofsson (1999); M98 Monet et al. (1998); N86 Nyman et al. (1986); 098 Olofsson et al. 
(1998); S95 Stanek et al. (1995); Y95 Young (1995); vLOO van Langevelde et al. (2000). 


stars and proto-PNe (PPNe) on the IRAS colour-colour di¬ 
agram along a curve corresponding to increasingly thicker 
and cooler CSEs. Jura (1993) suggest Mira periods continue 
to lengthen on the AGB but Whitelock et al. (1994) present 
evidence that Miras with P 400 days are a separate pop¬ 
ulation from longer-period objects. Semi-regular variables 
(SRa,b) are generally reported to have lower mass loss rates 
than Miras and could be stars entering (Kerschbaum & 
Hron 1992) or leaving (Young et al. 1993) the AGB. 

It is generally agreed that the stellar pulsations levitate 
the photosphere (Bowen 1988) until it is cool enough for 
dust grains to nucleate and grow at a few R t , as described by 
Gail & Sedlmayr (1998a, 1998b, 1999). Radiation pressure 
on dust then drives the wind away from the star via collisions 
with the gas. 

The CSEs of O-rich stars produce masers including 
SiO, H 2 0 22 -GHz, OH mainline and OH 1612-MHz emis¬ 
sion (Cohen 1989). High-resolution mapping of the winds 
from red supergiants (RSG) of M* ^ 10 Mq shows discrete 
H 2 0 vapour clouds 5 —15 au in radius which are 1 — 2 orders 
of magnitude denser than the surrounding wind (Richards 
et al. 1999). No direct measurements of the unbeamed sizes 
of H 2 O maser clouds around Miras have previously been 
published. The H 2 0 maser regions of AGB stars at distances 
of a few 100 pc are ~ 100 mas in radius. 

We observed four nearby AGB stars, IK Tau, U Ori, 
RT Vir and U Her, at 22 GHz using MERLIN. Its base¬ 
lines correspond to angular scales from about 10 to 200 


milli-arcseconds (mas), allowing us to image H 2 0 masers 
at high resolution and still detect all of the extended flux 
(Section 3.1). All four stars have 9.7-pm silicate features 
in emission (Olnon et al. 1986) and similar warm IRAS 
colours, log(F 25 /F 12 ) < 0, where F 25 and F 12 are the 25 and 
12 pm fluxes respectively. Other properties taken from the 
literature are given in Table 1. The absolute position of the 
star is given in column (2) and its velocity 14 with respect 
to the local standard of rest (Vlsr) is given in column (3). 
The period P is given in column (4). IK Tau, U Ori and U 
Her have regular periods and are classed as Mira variables; 
RT Vir is an SRb with no clear optical period. The optical 
magnitude range my is given in column (5) of Table 1. The 
distance d in column (6) is most reliable for RT Vir. The es¬ 
timates for IK Tau and U Ori using various methods agree 
with the results of the period-luminosity relationship within 
the uncertainties. For U Her this method gives double the 
distance derived from OH parallax. For all three Miras, the 
values of d are model dependent and for simplicity in all fu¬ 
ture calculations we use the Hipparcos distance of 133 pc for 
RT Vir and 266 pc for the other stars. The values quoted for 
the mass loss rates Mco in column (7) assume a gas-dust 
ratio of 120 — 200 and a fractional CO number density of 
(2 — 5) x 1CP 4 . The CO data also provide the maximum ex¬ 
pansion velocities and sizes of the CSEs in columns (8) and 
(9) respectively, see Section 4.3. In Section 2 we describe 
our observations and in Section 3 we present the observa¬ 
tional results. Our analysis of the kinematics and dynamics 
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Table 2. Details of MERLIN 22-GHz observations, including the 
restoring beams and noise levels in quiet channels for datacubes 
made at full resolution. 


Source 

(i) 

Date 

(yymmdd) 

(2) 

Integ 

time 

(hr) 

(3) 

Vlsr 
( km s _1 ) 

(4) 

©B 

(mas) 

(5) 

0 "rms 

(mjy 

brn" 1 ) 

(6) 

IK Tau 

940415 

11.0 

+34 

15 

10 

U Ori 

940417 

12.6 

-39 

15 

12 

RT Vir 

940416 

11.8 

+18 

20 

12 

U Her 

940413 

13.7 

-15 

15 

14 


of the circumstellar envelopes is given in Section 4. We in¬ 
vestigate the properties of individual water maser clouds in 
Section 5. We summarise our findings on maser clouds and 
the properties of AGB stellar winds in Section 6. 


2 OBSERVATIONS AND DATA REDUCTION. 

The four stars were observed in 1994 April using 5 tele¬ 
scopes of MERLIN: the Mk2 telescope at Jodrell Bank and 
outstation telescopes at Pickmere, Darnhall, Knockin and 
Cambridge. The maximum baseline was 217 km, giving a 
minimum fringe spacing of 13 mas at 22.235 GHz. Table 2 
gives further details of the observations, including the date 
and duration of the observations and the central velocity of 
the MERLIN spectral band Vlsr (columns (2), (3) and (4) 
respectively). Here and elsewhere radial velocities are given 
with respect to the Local Standard of Rest (LSR). We used 
a spectral bandwidth of 2 MHz divided into 256 frequency 
channels, which gave a velocity resolution of 0.1 km s -1 . 
The quasar 3C84 was observed for ~1 hr once or twice at 
the appropriate frequency for each source. It had a flux den¬ 
sity of 23 ± 1 Jy at that time (Terasranta, priv. comm.) and 
was used to calibrate the bandpass and set the flux scale for 
each source, with a final accuracy of ~ 10 percent. 

We reduced the data as described in Richards et al. 
(1999) for 22-GHz observations, using local MERLIN- 
specific programs and AIPS. LHC and RHC polarisation data 
were observed and calibrated separately and simultaneously, 
but all maps were made in total intensity. About 10 per¬ 
cent of the data were unusable, mainly due to bad weather. 
We derived corrections for phase and amplitude errors us¬ 
ing self-calibration only, so we could not obtain accurate 
absolute positions. We reweighted the calibrated data to at¬ 
tain the optimum combination of resolution and sensitivity 
and mapped and CLEANed each source using the AIPS task 
IMAGR. 

The FWHM (full width half maximum) of the restoring 
beam (0b) is given in column (5). The MERLIN beam at 
these declinations is moderately elliptical but using a circu¬ 
lar restoring beam of equivalent area did not introduce any 
artefacts and makes the maps easier to interpret. The typical 
rms noise for quiet channels, cr rms , is given in column (6). In 
the presence of bright emission (exceeding 1 Jy beam -1 ) the 
noise can be ~ 1 percent of the channel peak due to deconvo¬ 
lution errors arising from the sparse baseline coverage. Note 
that these parameters refer to the data cubes mapped at full 
velocity resolution (0.105 km s -1 ), which were used for all 


quantitative analysis of individual maser clouds. Additional 
maps were made at lower spatial and velocity resolution to 
study the overall properties of the H 2 O envelopes at greater 
sensitivity. Some of these maps are presented in Figs. 2 to 4. 

We fitted 2-D elliptical Gaussian components to each 
patch of emission in each channel in each datacube in order 
to measure the peak flux density 7 C , the position relative to 
the reference feature used for self-calibration, the total area 
and the total flux density S c . The relative position uncer¬ 
tainty <7p OS is proportional to the (dirty beam size) /(signal to 
noise ratio), as described analytically in Condon (1997) and 
Condon et al. (1998) and adapted for the MERLIN beam by 
Richards (1997) and Richards et al. (1999). This is typically 
1 mas for an isolated component with I c « 1 Jy bm _1 and 
0.1 mas for I c <; 10 Jy bm _1 . The FWHM of each compo¬ 
nent, s, was found by deconvolving the restoring beam from 
the component area, with an uncertainty cr s = \Z2a pos . 

The fitted components were grouped into features if 
three or more components with I c > 3<j rms occurred in ad¬ 
jacent channels with positions overlapping to within the po¬ 
sition error or component size. Non-matched components 
were discarded, as were any others which coincided with the 
beam sidelobe structure. 

3 OBSERVATIONAL RESULTS 

3.1 MERLIN spectra and maps 

Maps of the integrated 22-GHz emission from each source 
are shown in Figure 1. The emission comes in each case 
from elliptical regions ~0.3 arcsec in extent, in which are 
embedded hot spots ~20 mas in extent which contribute 
most of the flux. Table 3 lists for each source the total angu¬ 
lar extent of the emission detected, at the 5<r rms level A#tot 
(column(2)), together with the stellar velocity (column (3)), 
the minimum and maximum of Doppler velocity (columns 
(4) - (5)) and the range A Vlsr (column (6)), the peak flux 
density .Speak (column (7)) and the integrated emission Stot 
(column (8)). Figure 1 also shows estimated stellar positions 
derived in Section 4.1. It is clear that in each case the bright¬ 
est emission is mostly found towards the outer parts of the 
H 2 O envelope, with little emission from the direction of the 
star. 

Figs. 2 to 5 show the results in more detail: a MERLIN 
spectrum for each source, together with maps of the emis¬ 
sion integrated over velocity ranges of ~1.2 km s -1 . Individ¬ 
ual channel maps contain mostly well-separated, slightly re¬ 
solved patches of emission. Comparison with single dish data 
(autocorrelations and independent monitoring e.g. Rudnit- 
skij et al. 2000) shows that MERLIN detected all the flux 
to within the errors (< 10 percent). Note however that the 
extreme red-shifted emission from IK Tau extends right up 
to (and possibly beyond) the edge of our observing band. 
Henceforth red- and blue-shifted are used to denote emis¬ 
sion with Vlsr > V* and Vlsr < V* respectively. 

Fig. 6 shows the positions of the fitted maser compo¬ 
nents in each channel, with grey-scale used to indicate veloc¬ 
ity. In Section 5 we analyse the properties of the individual 
maser components. Here we concentrate on the global as¬ 
pects of the maser distributions around each star. Results 
for each source are given in the following four Sections. 
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Table 3. Global properties of the H 2 O maser envelopes. 


Source 

Af?tot 

V* 

Vlsr 

range 

AVlsr 

“Speak 

Stot 


(mas) 

(km s -1 ) 

min 

max 

(km s -1 ) 

(Jy) 

(Jy km s -1 ) 




(km 

s- 1 ) 




(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

IK Tau 

445 

+34.0 

+21.2 

> +46.7 

> 25.5 

25 

>108 

U Ori 

235 

-39.5 

-41.5 

-34.8 

6.7 

26 

20 

RT Vir 

370 

+18.2 

+9.5 

+26.6 

17.1 

394 

755 

U Her 

300 

-14.5 

-22.4 

-9.8 

12.6 

22 

45 




200 150 100 50 0 -50 -100 

Right Ascension Offset (mas) 



200 150 100 50 0 -50 -100 
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Figure 1. Integrated H 2 O maser emission from four AGB stars observed by MERLIN in April 1994. The origin is at the position of the 
reference feature. The cross shows the estimated stellar position (Section 4.1). The restoring beam is shown in the lower left corner. The 
grey-scale shows all emission, the contour plots were made after blanking each channel at its 5cr rms level. The contour levels are at (1, 
2, 4...)x the lowest contour level. This is at 0.1, 0.05, 0.3 and 0.3 Jy km s -1 beam -1 for IK Tau, U Ori, RT Vir and U Her respectively. 
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Figure 2. Top MERLIN spectrum of IK Tau 22-GHz emission. Bottom Maps of emission integrated every 12 chan 
the position of the reference feature. Each pane is labelled with the central velocity. The contour levels are (1, 2, 4 
The restoring beam is shown in the first pane. The panes containing the brightest emission (at 46.4, 45.1 and 31.2 k 
artefacts due to dynamic range limitations. 
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Figure 3. Top MERLIN spectrum of U Ori 22-GHz emission. Bottom Maps of emission integrated every 10 channels. The origin is at 
the position of the reference feature. Each pane is labelled with the central velocity. The contour levels are (1, 2, 4...) X 0.04 Jy bin -1 . 
The restoring beam is shown in the first pane. The pane containing the brightest emission (at -37.4 km s —1 ) contains some artefacts due 
to dynamic range limitations. 


3.2 IK Tau 

IK Tau lias the most complex spectrum and the widest ve¬ 
locity range of the four sources we studied. There are two 
bright spectral peaks asymmetrically distributed on either 
side of the stellar velocity: the blue-sliifted peak by 6 km s _1 
and the red-shifted peak by 12 km s -1 . The channel maps 
of IK Tau show that these bright emission peaks are due 
to multiple bright clumps spread over an elliptical region 
~ 200 mas E-W and ~ 100 mas N-S. Fainter clumps of 
emission are found in all channels, from an almost spherical 
region. Moderately blue-shifted emission to the west tends 
to be brighter than that to the east; the opposite is true 
for red-shifted emission. At velocities close to V* a clump 
to the south dominates the emission. The extreme red- and 
blue-shifted emission features are offset towards the east. 

Previous MERLIN observations by Yates & Cohen 
(1994) in 1985 detected only the bright elliptical region of 
emission but not the fainter spherical envelope. No individ¬ 
ual maser features from the earlier MERLIN maps can be 
matched to features in the new maps. On the other hand, 
Yates & Cohen (1994) found a close correspondence between 
MERLIN and VLA maps made 16 months apart (Lane et al. 
1987; Bowers et al. 1993). Many maser features had survived 


over that time period and Yates & Cohen (1994) were able 
to measure expansion of the H2O maser envelope. 

Although the overall appearance of the envelope in 1994 
is of a spherically symmetric shell, nevertheless there is a 
clear offset between the near and far sides of the shell, 
with red-shifted emission mainly to the east and blue-shifted 
emission mainly to the west. This is particularly noticeable 
in Figure 6. The E-W velocity segregation can be explained 
if the shell has an equatorial density enhancement so that 
the brightest masers lie in an oblate spheroid. The plane 
of the equator is at an angle of inclination i to the line of 
sight, with the eastern end of the polar axis approaching us. 
Hale et al. (1997) report that 3.1 (im speckle interferometry 
shows an approximately circular disc of warm dust 140 mas 
in diameter with a resolution of 40 mas. If this is the pro¬ 
jection of the inner regions of the H2O maser shell, with a 
minimum axial ratio within the uncertainty, this constrains 
i 45°. Bowers (1991) fig. 2 shows model maps for an oblate 
spheroid at * = 45° producing an asymmetric appearance 
at velocities around midway between V, and the extreme 
velocities. This model does not include acceleration, which 
produces brighter tangential beaming. Kirrane (1987) found 
that the OH mainline masers trace a biconical outflow with 
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Figure 4. Top MERLIN spectrum of RT Vir 22-GHz emission. Bottom Maps of emission integrated every 12 channels. The origin is at 
the position of the reference feature. Each pane is labelled with the central velocity. The contour levels are (1, 2, 4...) X 0.07 Jy bm _1 . 
The restoring beam is shown in the first pane. The panes containing the brightest emission (at 23.2, 13.1 and 11.8 km s —1 ) contain some 
artefacts due to dynamic range limitations. 


a wide opening angle, surrounded by OH-1612 MHz masers 3.3 U Ori 
in the equatorial region, with i « 60°. 


U Ori has the narrowest velocity range of the four sources, 
with a single dominant spectral peak at a red-shift of 
2 km s -1 and mainly red-shifted emission. The brightest 
emission, at —37.4 km s -1 , comes from an extended region 
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l n 



Figure 5. Top MERLIN spectrum of U Her 22-GHz emission. Bottom Maps of emission integrated every 10 channels. The origin is at 
the position of the reference feature. Each pane is labelled with the central velocity. The contour levels are (1, 2, 4...) X 0.07 Jy bin -1 . 
The restoring beam is shown in the first pane. The panes containing the brightest emission (at —14.6, —15.6 and -16.7 km s” 1 ) contain 
some artefacts, mostly near the centre, due to dynamic range limitations. 


to the south west, while the remainder of the emission is to 
the north and east. The distribution as a whole is elongated 
NE-SW, at a position angle of ~30°. The maser envelope is 
the most sparsely filled of the four we studied (Figure 6). 

U Ori has been monitored at 22 GHz by Rudnitskij 
et al. (2000) since 1980 using the Pushchino single dish ra¬ 
dio telescope. They observed H2O maser flares, which they 
propose arise due to shocks propagating radially out from 
the star. Spectacular OH maser flares have also been ob¬ 
served (Pataki & Kolena 1974; Chapman & Cohen 1985 and 
references therein). 

Bowers & Johnston (1994) mapped U Ori in 1988 using 
the VLA and found a similar elongated distribution at po¬ 
sition angle ~60°. The major difference between our maps 
and theirs is that they detected weak blue-shifted emission 


out to -49 km s -1 near the centre of the distribution, close 
to our estimated stellar position. The position angle of the 
H2O envelope observed by Bowers & Johnston (1994) in 
1988 is identical to the position angle of the OH 1612-MHz 
maser flare measured by Chapman & Cohen (1985) in 1983. 

3.4 RT Vir 

RT Vir has the greatest 22-GHz peak flux density and to¬ 
tal flux of the four sources, with two bright spectral peaks 
~5 km s' 1 each side of the stellar velocity. 

The emission has an angular extent < 370 mas, con¬ 
taining an inner ring of bright masers ~ 150 mas in diame¬ 
ter. As in IK Tau the emission contributing to the strongest 
spectral peaks is due to many compact maser regions spread 
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Figure 6. The positions of the H 2 O maser components for each source. Each panel is labelled with the source name. The origin is at 
the estimated stellar position. The diameter of each spot is proportional to the logarithm of its flux density. The paler spots correspond 
to more blue-shifted emission and the darker spots correspond to more red-shifted emission. 


over a large region. The total extent of the emission region is 
greatest near the stellar velocity and least near the extreme 
velocities, as expected for an expanding spherical shell, but 
again as in IK Tau the extreme red- and blue-sliifted fea¬ 
tures are displaced from each other. Extreme blue-shifted 
emission from RT Vir comes from near the centre of the 


channel maps, but the extreme red-shifted emission is to 
the south east. The most extended emission consists of faint 
patches to the N (slightly blue-shifted) and to the south 
(slightly red- and blue-shifted). 
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The integrated emission shown in Figure 1 gives the 
impression of a well-filled spherically symmetrical shell, but 
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in fact there is a systematic displacement between the near 
and far sides of the envelope, with most red-shifted masers 
to the east and blue-shifted masers to the west (Figure 6 
This E-W velocity segregation can be explained if the shell 
has an equatorial density enhancement so that the brightest 
masers lie in an oblate spheroid. 

Previous maps of RT Vir show a similar displacement 
between the near and far sides of the H2O envelope (Imai 
et al. 1997, Bowers & Johnston 1994, Yates & Cohen 1994, 
Bowers, Claussen & Johnston 1993). Despite this similarity, 
individual maser features cannot be matched. 

RT Vir has been monitored at 22 GHz at Pushchino 
since 1985 (Mendoza-Torres et al. 1997; Lekht. et al. 1999). 
It is extremely rapidly variable, with individual flares lasting 
less than 3 months. The MERLIN spectrum has a peak 150 
Jy brighter than measured at Pushchino a month earlier. 

RT Vir appears to have an equatorial density enhance¬ 
ment in its 22-GHz maser distribution such that the bright¬ 
est emission comes from an oblate spheroid tilted in the 
plane of the sky so that the western side is approaching us. 
This structure persists much longer than the lifetime of iden¬ 
tifiable individual maser clouds. The model maps of Bowers 
(1991) fig. 2 qualitatively illustrate the appearance of the 
brightest emission (as for IK Tau). The polar axis is likely 
to be at a projected angle similar to the direction of the 
red-blue velocity offset. 

3.5 U Her 

The emission from U Her is strongest near the stellar ve¬ 
locity and mainly blue-shifted. The emission comes from an 
elliptical ring-like region centred on our estimated stellar 
position, and elongated roughly north-south. The brightest 
emission lies to the south and west. The emission from each 
channel map is likewise ring-like, with no systematic change 
of radius with velocity. The most complete rings are seen 
near the stellar velocity, whereas at the extreme velocities 
just a single emission region is seen, offset from the stel¬ 
lar position by ~80 mas in each case. The most red-shifted 
emission occurs near the inner edge of the ring towards the 
west. The extreme blue-shifted emission lies to the south 
south west. The only emission clearly outside the ring is a 
moderately blue-shifted spur ~ 100 mas to the west. 

A similar ring-like distribution was observed by Bowers 
& Johnston (1994) in 1988 with the VLA, when the source 
was an order of magnitude brighter. However there is no de¬ 
tailed correspondence of individual maser features between 
the two epochs. The extreme blue-shifted emission which 
formed the north-eastern part of the ring in 1988 has now 
faded. The radius of the ring may appear to have increased 
by ~20 mas in the 5 years between the two set of observa¬ 
tions, but this result is only marginally significant given the 
70-mas beam of the VLA and the lack of detailed agreement 
at the level of individual components. Colomer et al. (2000) 
also found a ring distribution using the VLA in 1990, but 
again with no detailed correspondence between individual 
maser features at the different epochs. 

Previous MERLIN measurements by Yates & Cohen 
(1994) in 1985 detected only a single extended region of 
emission near the stellar velocity. 

The strong differences in appearance at different epochs 
show that individual clouds in the envelope of U Her cannot 


Table 4. Estimated stellar positions ( X , Y) with respect to the 
reference feature for each source. 


Source 

(A, V) 
(mas) 

Error in X and Y 
(mas) 

IK Tau 

(-38, 67) 

4 

U Ori 

(IT 57) 

9 

RT Vir 

(32, -40) 

3 

U Her 

(65, 7) 

5 


be followed for more than ~2 stellar periods. This is consis¬ 
tent with the shock model discussed later in Section 5.4. 
U Her, like U Ori, has also undergone OH maser flares 
(Etoka & Le Squeren 1997), which may be another result 
of shocks. 


4 THE DYNAMICS OF THE 

CIRCUMSTELLAR ENVELOPES 

4.1 Stellar position 

Much of our further analysis requires an estimate of the 
stellar position. We used two methods to estimate the centre 
of expansion ( X , Y ), assumed to be the stellar position. 
Both methods assume that V* is accurately known and use 
least-squares minimisation routines. 

The shell-fitting method of Yates (1993) Ends the point 
which is most nearly equidistant from all components in an¬ 
gular separation and in velocity. This assumes the maser 
components have a point-symmetric distribution about the 
centre of a spherical shell. However for non-spherical shells 
it is more reliable to use a method based on the quench¬ 
ing radius (Section 5.3). Values of (A', Y) are found which 
maximise the angular separation from all components with 
velocities within the range V* ± (AVlsr/u), where n is be¬ 
tween 4 and 8 (see Table 3). 

In practice, similar results were obtained using both 
methods for IK Tau, RT Vir and U Her. For U Ori the 
shell-fitting method did not converge, but the quenching ra¬ 
dius method gave consistent results for different values of n. 
Table 4 gives (A, Y) found by the quenching radius method. 
The error quoted is the difference between this position and 
the mean position found by all methods. 

Figs. 2 to 5 ( bottom ) show that in all four CSEs either 
or both of the extreme blue- and red-shifted features are 
offset from the assumed stellar position by 20 115 mas. 

We did not consider the mid-point of the extreme blue- and 
red-sliifted features as the stellar position as in every case 
this would be very asymmetric with respect to most of the 
masers. The offsets could be explained by turbulence of 1.1 
- 1.6 km s" 1 (Chapman et al. 1994; Chapman 1985), by 
dumpiness, or by more systematic causes discussed in Sec¬ 
tions 3.1 and 5.4. 

4.2 Kinematics of the H 2 O shells 

Previous MERLIN observations of circumstellar H2O masers 
have suggested a thick-shell model, in which the masers 
are distributed irregularly but the velocity Held is regu¬ 
lar and spherically symmetric, increasing steadily with ra- 
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Figure 7. Radius-velocity plots of the 22-GHz masers around IK Tau, U Ori, RT Vir and U Her, showing the angular separation a from 
the assumed stellar position, as a function of Vlsr- The diameter of each spot shown proportional to the logarithm of its flux density. 
The dashed lines show ellipses centred on V - * which were fitted by eye to the inner and outer limits of the maser distribution. The dotted 
line shows an alternative inner limit for IK Tau (see Section 5.3). 


dius through the H 2 O maser region (Yates & Cohen 1994, 
Richards et al. 1999 and references therein). That model also 
describes the present observations well. In Fig. 7 we show 
the radius-velocity plots for the four sources. The inner and 
outer edges of the H 2 O envelope were estimated by fitting 


the dashed ellipses by eye to the data. From these we obtain 
the inner and outer radii of the H 2 O maser zone, which we 
denote by n and r 0 , together with the corresponding ex¬ 
pansion velocities Vi and v 0 . These parameters are listed in 
columns (2) - (7) of Table 5. The error in the estimated 
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stellar position (Section 4.1) is applicable to the estimates 
of n and r 0 and is given in column (8) of the Table. 

Two alternative inner limits are shown for IK Tau, since 
the limit (a) with the larger value of ?'i is well defined but a 
few features lie within this giving a smaller value (b). The 
limits for U Ori and U Her are more uncertain but have been 
informed by the appearance of these stars at other epochs 
(e.g. Bowers & Johnston (1994), Murakawa priv. comm.) 
when emission appeared at different position angles and/or 
velocities. 

In all cases the expansion velocity is less than the escape 
velocity at the inner edge of the H 2 O maser zone, but easily 
exceeds the escape velocity at the outer edge of the H 2 O 
maser envelope. This was already established by Yates & 
Cohen (1994) for IK Tau and RT Vir, and seems to be a 
characteristic of circumstellar H 2 O envelopes in general. 

The shell limits were used to derive the logarithmic ve¬ 
locity gradient e = d(ln v)/d(lnr) in each circumstellar en¬ 
velope. These are listed in column (9) of Table 5, along with 
the error cr e , the linear velocity gradient K v and its uncer¬ 
tainty <tk v (in columns (10) to (12) respectively). These are 
average values for the H 2 O maser shells as there are local 
variations in each CSE. Acceleration is strongest nearer the 
star in each case. 

The values of e are similar in the four stars: they range 
from 0.5 to 0.9, and are very similar to the values derived 
from MERLIN observations of supergiants (Yates & Cohen 
1994 and references therein). These are the first direct mea¬ 
surements of t for low-mass AGB stars, since ours are the 
first measurements to resolve the thickness of the H 2 O shells. 
Previous MERLIN measurements by Yates & Cohen (1994) 
made at lower angular resolution (without the Cambridge 
telescope) could only set upper limits on t for IK Tau and 
RT Vir. 

The total H 2 O maser luminosity 4? (equivalent isotropic 
luminosity) is given in column (13) of Table 5 for complete¬ 
ness. Although the three Miras all have peaks of similar 
intensity (Table 3) and are at a similar distance, U Ori has 
only about a quarter of the luminosity of IK Tau, U Her 
being at an intermediate value. RT Vir is almost twice as 
luminous as the brightest of the Miras. 

The present measurements are sufficiently detailed to 
allow us to investigate the variation of maser emissivity with 
radius in the CSEs. Figs. 7 and the spectra in Figs 2 - UHER 
show that more and brighter masers were detected in the 
inner parts of the shells. This is broadly consistent with 
the models of Cooke & Elitzur (1985) who predict a sharp 
increase in maser brightness near n and a gradual decline 
towards r 0 . 

The 22-GHz maser shells around IK Tau and RT Vir are 
well-filled enough for us to estimate the photon luminosity 
as a function of distance from the star. We used the model of 
a spherically symmetric velocity field within the shell limits 
given in Table 5 to assign each maser component a full set of 
vectors to describe its position and velocity in all 3 spatial 
dimensions, as described in Murakawa et al. (2002, MNRAS, 
submitted). We then calculated the flux density in sub-shells 
with a cross-section in the plane of the sky of 10 mas, and 
used this to find the 22-GHz photon rate per unit radius 
d&/dr. This is shown in Fig. 8. The results are robust at 
the 10 percent level to departures from spherical symmetry 
by up to a factor 2 in t (Richards et al. 1999). 



Figure 8 . The 22-GHz photon rate per unit radius for IK Tau 
and RT Vir, as a function of radial distance from the assumed 
stellar position. The bold numbers above the abscissa are esti¬ 
mates of the time taken for material in the stellar wind to reach 
that distance, see Section 5.3. 


The results shown in Fig. 8 were compared with the 
model results of Cooke & Elitzur (1985). Initially we used 
the mass loss rates A/co given in Table 1. For IK Tau d.^/dr 
reaches a maximum value of ~ 2 x 10 33 photons s -1 m -1 
at a radius of 6 x 10 12 m. This distance is within the range 
predicted by the models of Cooke & Elitzur (1985) for the 
given A/co- However, for RT Vir, the maximum d<fr/dr ~ 
4 x 10 33 photons s -1 m” 1 occurs at a radius of 3 x 10 12 m, 
which is ~ 3 times the radius predicted by Cooke & Elitzur 
(1985) for the given A/co- 

We return to this comparison later in Section 5.4, when 
we discuss the over-density of the H 2 O maser clouds with 
respect to their surroundings. 


4.3 The velocity fields of the stellar winds 

Fig. 9 shows plots of expansion velocity vs. radial distance 
for the four stars, compiled using our H 2 O data together 
with data on other molecular species taken from the litera¬ 
ture. All four stars have bright OH mainline masers which 
have been mapped (Chapman et al. 1991; Chapman et al. 
1994; Bowers & Johnston 1988; Richards et al. 2000; Kir- 
rane 1987), while IK Tau was mapped at 1612 MHz by Kir- 
rane (1987). The OH 1612-MHz maser flares in U Ori were 
imaged by Chapman & Cohen (1985) but as there was no 
clear shell structure there was no obvious way to assign a 
radial distance to the emission. CO data were taken from 
the references in Table 1. Only RT Vir is well resolved in 
CO (Stanek et al. 1995). The values given for the Miras are 
based on different models used by the various authors, but 
they also observed some of the same stars as Stanek et al. 
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Table 5. Parameters of the H 2 O maser envelopes (see Sections 4.2 and 4.3 for details). 


Source 

T*i 


V[ 

r 0 


Vo 

err 

e 

Ce 

A' v 

°K V 

4> 


(mas) 

(au) 

(km s -1 ) 

(mas) 

(au) 

(km s -1 ) 

(au) 



(10“ 12 


(10 42 











km s —1 

m -1 ) 

photons s _1 ) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(U) 

(12) 

(13) 

IK Tau(a) 

60 

16.0 

4.5 

250 

66.5 

16.0 

1.0 

0.89 

0.08 

1.52 

0.06 

4.63 

IK Tau(b) 

25 

6.7 

5.5 

250 

66.5 

16.0 

1.0 

0.46 

0.09 

1.17 

0.04 

4.63 

U Ori 

37 

9.8 

2.5 

120 

32.5 

7.5 

2.4 

0.92 

0.24 

1.48 

0.23 

0.85 

RT Vir 

45 

6.0 

4.0 

185 

25.6 

10.0 

0.4 

0.65 

0.05 

2.15 

0.08 

7.88 

U Her 

48 

12.8 

5.0 

175 

46.6 

9.5 

1.1 

0.50 

0.04 

0.89 

0.05 

1.92 


(1995). Comparison of the results for objects in common 
suggests a CO diameter of 10000 au for IK Tau, and sizes a 
factor of 3 — 5x smaller than those given for U Her and U 
Ori. The error bars for CO data in Fig. 9 reflect this. The 
distances were taken as 133 pc for RT Vir and 266 pc for 
the other stars. 

Fig. 9 show several striking features: 

• There is a general increase of expansion velocity with 
distance from the star, with some scatter. The strongest 
acceleration is in the region out to 30 au in the smallest 
CSE, RT Vir, and 70 au in the largest, IK Tau. 

• IK Tau and U Her also show evidence for gentler accel¬ 
eration continuing further outwards to ~500 au, and perhaps 
as far as the CO region at ~3000 au. 

• The best filled 22-GHz maser shells, around IK Tau and 
RT Vir, attain the greatest H 2 O velocities and contain faster 
emission at a given radius than do the shells of U Ori and 
U Her. 

• In U Ori, U Her and RT Vir, the OH mainline masers 
overlap the H 2 O shells in angular separations and velocities. 
However the H 2 O masers appear to have greater velocities 
than OH masers at a similar radius. 

Chapman & Cohen (1986) observed acceleration of the 
wind from VX Sgr out to at least 100 R * and showed that 
could be explained by an increase in the dust absorption ef¬ 
ficiency kd during the outflow. A similar kinematic pattern 
has been observed in other RSG e.g. S Per (Richards et al. 
1999) but hitherto it has been hard to confirm similar be¬ 
haviour in the smaller CSEs of Miras. Fig. 9 and Table 5 
show that acceleration occurs out to many tens of R* in 
these stars, and much further for two of the stars: IK Tau 
and U Her. 

If dust properties are assumed to be constant Elitzur 
& Ivezic (2001) predict that some acceleration does occur 
due to optical depth effects but terminal velocity is attained 
at ^ 10.R* and wind speed increases monotonically with M 
for silicate dust (Netzer & Elitzur 1993). Monnier, Geballe 
& Danchi (1998, 1999) found periodic and long-term fluc¬ 
tuations in the shape and intensity of the 9.7/im silicate 
feature towards IK Tau, U Ori and U Her. This suggests 
changes in the dust optical constants; possible mechanisms 
include annealing or solid diffusion (Gail & Sedlmayr (1999) 
and references therin), which could improve the absorption 
efficiency at greater r. 

Multi-epoch monitoring is continuing to investigate 
the effects of stellar variability on the velocity field in 
CSEs. Engels et al. (1997) found strong variability (“mode¬ 
switching”) in the 22-GHz line shape and velocity width 


towards OH 39.7+1.5, with no evidence for velocity gradi¬ 
ents at large radial distances. However, unlike the stars in 
our sample, this object has cool IRAS colours and a 9.7- 
/rm silicate feature in absorption. Engels et al. (1988) and 
Szymczak & Engels (1995) also found that H 2 O masers at¬ 
tain terminal velocity in objects with high M, whereas H 2 O 
masers in thinner-shelled CSEs do not. This is consistent 
with our results and suggests that stronger acceleration is 
favoured when starlight can more easily penetrate the CSE. 

5 THE PROPERTIES OF MASER CLOUDS 

5.1 Maser component sizes and temperatures 

We now analyse the properties of individual water vapour 
maser clouds. We begin by considering the Gaussian compo¬ 
nents fitted to each channel map. In the following sections 
we will deal with the maser clouds themselves. Fig. 6 shows 
the positions of the fitted components in each channel rela¬ 
tive to the estimated stellar position (Section 4.1) for each 
source. A complete list of the fitted components and their 
fitted parameters is available electronically via CDS. 

In each source, the majority of components are resolved. 
Table 6 summarises component properties. The apparent 
size of an individual maser component s (measured as de¬ 
scribed in Section 2) is the FWHM of beamed maser emis¬ 
sion from molecules with velocities within the 0.105 km s _1 
channel width, and is smaller than the physical size of the 
emitting region. The greater the maser amplification factor, 
the smaller the beaming angle (Elitzur 1992), and so under 
comparable conditions brighter maser components appear 
smaller. Column (2) gives the total number N c of 22-GHz 
maser components in each source and column (3) gives the 
fraction of components F c for which s > a s , where a s is the 
uncertainty in s. The error-weighted mean and scatter of s, 
s and <Tg are given in columns (4) and (5). The scatter is 
intrinsic to the data, not simply due to measurement un¬ 
certainties. Note that if RT Vir is at half the distance of 
the other sources, s represents a similar actual size in all 
sources apart from U Ori, whose components appear signif¬ 
icantly larger. Imai et al. (1997) resolved out well over half 
the flux from RT Vir using VLBI, but established that the 
brighter spots detected (mostly > 100 Jy in the total power 
spectrum) had a spatial FWHM of 0.3-1.2 mas, which is 
consistent with our results. 

The brightness temperature Tb c of each component was 
calculated from its total flux density S c and its size s. The 
minimum Tb c measured for each source was 10 5 K; the max- 
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Figure 9. Expansion velocity as a function of radial distance 
from the four stars IK Tau, U Ori, RT Vir and U Her. The H 2 O 
data are from Table 5. The OH data are from references given in 
Section 4.3. The CO values are outer limits from references given 
in Table 1. 
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Figure 11. Maser component FWHM s as a function of bright¬ 
ness temperature Tb c - The lines show error-weighted fits with 
slope cx as defined in Equation 1. 

ima are given in column (6) of Table 6. For a spherical mas- 
ing region, the relationship between beamed component size 
and brightness temperature can be approximated by 

socalogTec- (1) 

Elitzur (1992) develops a more complete model giving — 1 ^ 
a ^ —2 for spherical clouds, explaining how lower values 
correspond to more saturated maser emission if the other 
relevant factors are unchanged. 

Fig. 11 shows s as a function of logTec- We used error- 
weighted least-squares fits to find the slope a , which is given 
in column (7) of Table 6. For IK Tau, U Her and RT Vir, 
the results are broadly consistent with our simplified model 
and imply that the masers are saturated. The formal errors 
in a are < 0.2 and values of a < —2 may be due to the 
intrinsic scatter in s and to deviations from our assumptions. 
However, for U Ori a > 0 showing more drastic revision is 
needed, and this is discussed in Section 5.4. 

5.2 Measurements of maser clouds 


Figure 10. The labelled points show the measured values of r; as 
a function of (Mco) 2 ^ 3 /vi- The solid line is the predicted r q for a 
similar range of Mqq at a typical v\ = 5 km s —1 and / = 4x 10 —4 . 


Fig. 6 shows the grouping of components into necklace-like 
strings, with similar velocities and positions. Each string 
is thought to be produced by a single maser cloud, with 
the emission in each frequency channel picking out the line 
of greatest amplification within the velocity range of that 
channel. IK Tau contains a few clouds with distinct velocity- 
position gradients, but most have an apparently random dis- 
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Table 6. The error-weighted statistics of maser components and 
their brightness temperatures, see text for details. 


Source 

(i) 

Ac 

(2) 

F c 

(3) 

s 

(mas) 

(4) 

°5- 

(mas) 

(5) 

Max. T Bc 
(K) 

(6) 

a 

(7) 

IK Tau 

1490 

0.60 

3 

2 

5 x 10 9 

-4 

U Ori 

95 

0.82 

9 

3 

2 x 10 s 

+7 

RT Vir 

847 

0.73 

6 

2 

1 x 10 11 

-2 

U Her 

282 

0.55 

3 

2 

1 x 10 8 

-3 


tribution of components with velocity. In U Ori, only two 
clouds have any clear structure. However, in U Her much 

o 

of the blue-shifted emission is less random and in RT Vir g 
the majority of clouds have smooth (not necessarily linear) ° 
gradients, which are also seen in Fig 7. ° 

The measurements of each maser cloud derived from its jj> 
constituent components are given in Tables available elec- § 
tronically from CDS. Columns (1) and (2) give the flux- z 
weighted mean Vlsr of the cloud (Vlsr) and its total veloc¬ 
ity width (AC). Columns (3) and (4) respectively give the 
FWHM (AVi/ 2 ) of a Gaussian curve fitted to the cloud ve¬ 
locity profile and the uncertainty in this (ctavi/ 2 )- Columns 
(5)-(12) give the position of the error-weighted centroid of 
the cloud with respect to the assumed stellar position (Sec¬ 
tion 4.1): x, (jx , y and a y are the offsets and uncertainties 
in Cartesian coordinates and a, cr a , 6 and erg are in polar 
coordinates. Columns (13) and (14) respectively give the to¬ 
tal angular extent of the cloud (L) and its uncertainty (ox); 
column (15) gives the peak flux density (I), column (16) 
gives the total flux in the cloud (S') and column (17) gives 
the peak brightness temperature (7b). 

Using a distance of 133 pc for RT Vir and 266 pc for 
the other sources, cloud angular size L was converted to 
the cloud diameter l, which is likely to represent at least 
85 per cent of the true unbeamed size (Richards et al. 1999). 

Figs. 12 and 13 show the distribution of l, AV 1/2 and 
AV for each source. Individual cloud measurements of ac¬ 
curacy < 1.5 a (mostly for clouds comprising only 3 com¬ 
ponents) are not included. Table 7 gives the properties 
of an average cloud around each star derived from error- 
weighted means of the individual cloud measurements above 
the threshold of 1.5cr. Column (2) gives the total number N 
of clouds in the source. Column (3) gives the number of 
clouds N(l, 1.5) with l/ai > 1.5, where a\ is the uncertainty 
in l. The mean (7) and dispersion (cry) of their diameters are 
given in columns (4) and (5) respectively. Columns (6) gives 
the number of clouds lV(AVi/2l.5) with AV 1/2 > 1.5frAv 1/2 - 
These were used to find the mean (AV 1 / 2 ) and dispersion 

(cr-^y—) of the velocity FWHM and of the total velocity 
1/2 _ 

width (AV) and the uncertainty in this (tr-^y-), given in 
columns (7) - (10). 

Table 7 shows that the clouds around IK Tau have a 
significantly larger 7 and a smaller AV than those around 
RT Vir. U Ori and U Her have intermediate values. The 
IK Tau clouds also have the smallest AVl/ 2 - Fig- 12 shows 
that IK Tau contains a greater number of larger clouds than 
the other sources. Fig. 13 shows that RT Vir and U Her 
have more clouds with AVi/ 2 > 0.7 km s -1 than the other 
sources although RT Vir has almost as many with 0.4 < 



J_I_I_,_,_ 1 

5 10 


Cloud Diameter l (au) 

Figure 12. The distribution of diameters l of maser clouds in 
each source, for clouds with l > 1.5cri. 


AV 1/2 < 0.5 km s -1 , which is the peak range for IK Tau. 
The peak of the distribution of AV is at 1.6 — 2.0 km s -1 for 
RT Vir, three times the value for the other sources. RT Vir 
and IK Tau have a few clouds with AV ss 3AVth where AVth 
is the thermal line-width of H 2 O at 1000 K, ~ 1.4 km s -1 . 

The necklace-like series of components in Fig. 6 appear 
to be in random directions. We searched for systematic ve¬ 
locity gradients within maser clouds as follows. For each 
cloud we attempted to fit a straight line to u = Vlsr— V* as a 
function of position to find the gradient u! , taking the good¬ 
ness of fit as well as measurement errors into consideration 
in finding the uncertainty ct u /. Columns (11), (12) and (13) 
of Table 7 give the number of clouds N(u') with u' > 2a u i 
used to find the mean (it') and the maximum (Max. it') of 
it'. The minimum measurable it' was 3 x 10 -12 km s -1 m -1 . 
The typical velocity gradient v! was ~20 x 10~ 12 km s _1 
m -1 . This is at least an order of magnitude larger than the 
systematic velocity gradient K v (Table 5) within each ex¬ 
panding H 2 0 envelope. The most likely cause of this is local 
turbulence, which could be thermal, as AV is of the same 
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Figure 13. The distribution of the total velocity extent AV (dot¬ 
ted lines) and FWHM AV) / 2 (solid lines) of maser clouds in each 
source for clouds with AV 1/2 > 1-5ctav 1 / 2 ' 


order as AVth, or due to shocks or differential acceleration 
between regions of different density or dust:gas ratio (see 
Section 5.3). 

RT Vir has the strongest velocity gradients, with u' 
more than double the u' seen for the other stars. 

We also decomposed u' into dtt/da and du/d$, but 
found no statistically significant preferred direction for the 
individual cloud velocity gradients. 

Fig. 14 shows that the distribution of peak Tb for clouds 
around RT Vir reaches a maximum in a range 10 x hot¬ 
ter than for the other stars, as well as extending to values 
lOOx higher. For unsaturated masers, AVi/ 2 narrows as the 
maser gets brighter and then rebroadens as the maser sat¬ 
urates (Elitzur 1992). Although Figs. 13 and 14 suggest a 
correlation between greater AV 1/2 and Tb, there is no clear 
relationship within the cloud measurements for each source. 



Peak Brightness Temperature T B (K) 


Figure 14. The distribution of the peak brightness temperatures 
of maser clouds in each source. 
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Table 7. The average properties of clouds around each source. See Section 5.2 for details. 


Source 

(i) 

N 

(2) 

JV( 11.5) 

(3) 

l 

( 

(4) 

°T 

au) 

(5) 

N{ AV 1/2 1.5) 

(6) 

AVi/a 

(km 

(7) 

ff AV 1/2 

S- 1 ) 

(8) 

(km 

(9) 

ff 4V 

s- 1 ) 

(10) 

N(u') 

(11) 

u' 

(10- 12 

(12) 

Max. u 

km s —1 m —1 
(13) 

IK Tau 

256 

49 

3.8 

0.2 

120 

0.46 

0.02 

0.84 

0.06 

31 

11 

68 

U Ori 

14 

10 

3.0 

0.8 

7 

0.68 

0.12 

1.07 

0.24 

3 

17 

45 

RT Vir 

55 

26 

2.2 

0.2 

52 

0.67 

0.04 

1.67 

0.11 

18 

43 

190 

U Her 

34 

24 

2.6 

0.3 

24 

0.78 

0.05 

1.03 

0.12 

9 

11 

23 


5.3 Cloud densities and survival 

The inner radius r; of 22-GHz H 2 O maser emission is de¬ 
termined by the quenching number density n q , above which 
the collision rate exceeds the radiative decay rate of the up¬ 
per level of the maser transition. At 1200 K, for a fractional 
H 2 O number density / = 4 x 10 _4 tih 2 (where riH 2 is the H 2 
number density), n q = 5x 10 15 (Yates et al. 1997). Cooke 
& Elitzur (1985) predicted that n oc M 2 ^ 3 fv\. This relation 
is plotted in Fig. 10, using the relevant values given in Ta¬ 
bles 1 and 5. Fig. 10 shows that the observed n of the four 
sources lie well above the inner radii inferred from A/co- 
This suggests the maser clouds are denser than average for 
a homogenous wind, as if they were produced by localised 
episodes of mass loss at a higher rate. 

The required mass loss rates (Mh 2 o, uncertainty 
cr M H 2 o^ are derived as described in Richards et al. (1999). 
These are given in column (2) and (3) of Table 8. The CO 
mass loss rates A/co are given in Table 1 and the ratio of 
the cloud density to the average density at a given radius 
is Mh 2 o/Mco = f P , given in column (4) of Table 8. The 
quenching density also implies that all maser clouds at a 
given temperature have the same density, regardless of the 
mass loss rate. 7 is used to find the volume filling factor for 
the clouds /v given in column (5) of Table 8. 

We adopt r p = n + (r 0 — n)/4 as a typical radius for 
the best-filled and brighest part of the shell, where the true 
cloud diameter is probably closest to the average measured 
cloud size 7. Following Richards et al. (1999), we estimated 
the cloud density at ?' p by extrapolation from n q and used 
this to find the typical mass of a cloud of diameter 7, A/h 2 o, 
given in column (6). 

The cloud sound crossing time t cs is given for 7 in col¬ 
umn (7). The time t s taken for a cloud to travel from n to 
r 0 is given in column (8) and used to find the apparent num¬ 
ber of clouds formed per stellar period Np, given in column 
(9). Fig. 8 shows the time since the clouds left the star as 
a function of radius, where the time at r 0 is time taken to 
reach n (assuming a constant velocity of v; from R * to n) 
plus t s . IK Tau appears to form more clouds per period than 
the other stars although the value for RT Vir is uncertain 
due to its irregular period. 

With tcs of about 1.5-3 yr, all of the sources have 
tcs ~ l/10fs making them vulnerable to turbulent disrup¬ 
tion. Indeed, observations at other epochs (see references in 
Sections 3.2 - 3.5) show that individual clouds can be iden¬ 
tified at intervals of up to 18 months but not after 2-5 yr, 
depending on the source. There is only indirect evidence at 
present as to whether this indicates that the clouds dissipate, 
or whether their maser emission in our direction simply de¬ 
creases. The small volume of the shell occupied by clouds 


(fv < 0.01), see Table 8, together with the long transit 
time for gas through the H 2 O maser zone compared with 
the sound-crossing time for a cloud (t s and t cs typically 20 
yr and 2 yr respectively) suggest that a change in the maser 
output towards us is the most likely explanation. Such a 
change could easily come about due to the large internal 
velocity gradients u' that we observe in random directions 
within individual maser clouds (Table 7 and Section 5.2). A 
direct proof of this hypothesis would be to observe clouds 
reappear at a later date. 

The ratio 7/r p is remarkably similar for all four stars, 
between 0.12 - 0.20, where the scatter is as expected from 
the measurement uncertainties. If the clouds originate from 
the star their birth radius is a fraction of the stellar radius 
given by 7(/7*/r p )/2/7* = 0.08 ± 0.02/7*. 

Returning to Table 8, the total mass in the clouds we 
observe in each CSE is given in column (10) by Mtot h 2 o = 
Mh 2 oN (N is given in Table 7). The total mass injected at 
Ti during t B is A/ to t inj = t s Mco, given in column (11). The 
uncertainty in quantities derived from our measurements are 
^ 50 per cent, and comparisons with A/co may be less ac¬ 
curate. 

Even allowing for such uncertainties, some trends ob¬ 
served in all four CSE are strongly significant. 

• The maser clouds are one or two orders of magnitude 
denser than the surrounding gas, showing that in all these 
stars, mass is lost in the form of density-bounded discrete 
clumps. 

• A/tot inj/A/tot h 2 o is between 1 and 5 (to within our ac¬ 
curacy of 50 per cent), consistent with the filling factor and 
overdensity of the H 2 O maser clumps. 

• For each star the radius of the clouds, extrapolated back 
to the stellar surface, is ~ 0.1/7*. 

Prior to this work, the size of Mira clouds had not been 
measured. However, we (e.g. Richards et al. 1999; Murakawa 
et al. 2002, MNRAS, submitted) have previously shown that 
RSG H 2 O maser clouds are overdense by a similar amount, 
and are a similar relative size (~ 10 per cent) when extrap¬ 
olated back to /?*, although the absolute values of l, /7*, n 
and A/ are an order of magnitude greater for RSG. However, 
as befits their larger size, RSG cloud diameters are greater 
than the velocity resonance length and they survive for at 
least 9 years (Richards Yates & Cohen 1996, 1998). 

These results suggest that the cloud size is not primar¬ 
ily determined by local physical or molecular processes in 
the wind such as shocks or cooling, as these would operate 
on the same scales around RSG as around AGB stars. It 
could be that clumps somehow grow between the star and 
n, so the larger M and r q of RSG means the clouds reach 
a larger mass. Alternatively, the cloud size could be directly 
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determined at the stellar surface, related to convection cells 
(possibly producing local chemical enrichment), pulsation ir¬ 
regularities producing ejections, or cooling above starspots 
(Frank 1995). 

5.4 Maser amplification 

Maser amplification occurs over a velocity resonance length 
51 where Ulsr changes by AVth. Inspection of the individ¬ 
ual channel maps (Figs. 2-4) shows that there is only oc¬ 
casional blending, and the low filling factors imply there is 
< 1 per cent chance of extra amplification due to one cloud 
overlapping another along the line of sight unless their dis¬ 
tribution is very non-random. The actual clouds could be 
any shape, but for simplicity we assume they are spher¬ 
ical (however see Section 5.4). Using the method applied 
to the RSG S Per by Richards et al. (1999) and the shell 
parameters in Table 5, we estimate SI for each source at 
r p in the limiting cases of a cloud along the line of sight 
to the star (radial beaming) and in the plane of the sky 
with the star (tangential beaming). We find for IK Tau 
6 £ 51 ;> 5, for U Ori 51 « 6, for U Her 10 £ 51 )> 5 
and for RT Vir 4 <; 51 <; 3 au. The radially-beamed 51 
exceeds the tangentially-beamed value by a larger amount 
for smaller values of the logarithmic velocity gradient e. In 
all cases 51 > l (Table 7), so the velocity coherence depth 
is not limiting the amplification path, unless there is some 
systematic asymmetry in the distribution of cloud shapes. 

If the only force on clouds is due to radiation pressure on 
dust, their diameter in the tangential and radial directions 
will increase proportional to r/r\ and (r/n) 6 respectively. 
If the clouds are on average initially spherical, they will re¬ 
tain their shape where t is not much less than 1. However, 
the clouds near r\ could be flattened by shocks due to the 
stellar pulsations. This favours tangential 22-GHz emission, 
perpendicular to the direction of the shock, not only due 
to the longer path length for maser amplifeation, but also 
because the pump IR photons escape along the short axis 
(Elitzur et al. 1992). In this geometry, the beaming angle is 
independent of Tb and the observed component angular size 
is equivalent to the projected width of the flattened cloud 
(Elitzur 1992). 

Thus, shock compression would explain the large scale 
of emission inferred from the visibilities in the brightest 
channel of U Ori (Section 3.3) and its large values of a and 
s in Table 6. The ring-like appearance of the masers around 
U Her is also consistent with shocked clouds. For unsatu¬ 
rated masers, Tb is roughly exponentially proportional to 
a linear function of the path length. The brightest masers 
around U Her have Tb « 10 8 K. If this arises from clouds 
with an axial ratio of 2:1 viewed along the long axis, clouds 
viewed along the short axis will have Tb « 10 4 K, so clouds 
lying close to the line of sight to the star are well below 
the detection threshold. The average cloud diameter 7 was 
calculated on the basis of spherical clouds. If we only detect 
flattened clouds, l becomes the long axis, leading to an over¬ 
estimate of 7. However, if the clouds around the front and 
back caps are present but not detected, then N and Np are 
underestimated but the accuracy of other quantities given 
in Table 8 is not greatly affected. 

The 22-GHz photon luminosity rate from the whole of 
each shell, <f>, is shown in column (13) of Table 5. For the Mi¬ 


ras, the average d4>/dr is (0.3 — 0.6) x 10 30 photons s _1 m _1 
but for RT Vir it is 3 x 10 3 ° photons s _1 m _1 , i.e. ^ lOx 
higher. The individual masers around RT Vir also reach 
the highest Tb, although this object has the lowest M, the 
smallest shell and the smallest clouds. It also has the small¬ 
est, most irregular, stellar pulsations. It shows the largest 
velocity gradients ( u 1 and K v ) and it is the only source 
with AU > AVth (Table 7). The masers are very rapidly 
variable, suggesting highly efficient but unsaturated ampli¬ 
fication. Small clouds with a steep velocity gradient could 
allow the clouds to remain optically thin at the IR pump 
photon frequencies (Neufeld & Melnick 1991, Elitzur 1992) 
and achieve very efficient 22-GHz maser pumping. This is 
analogous to the suggested enhanced masing from shock- 
compressed clouds in U Her and U Ori. Shock compression 
alone is a less likely explanation in the case of RT Vir as 
its bright masers extend well beyond n and its optical pe¬ 
riod and amplitude are poorly defined and small (Table 1), 
suggesting the stellar pulsations are either weak or very com¬ 
plex. 

Another possiblity is that the wind from RT Vir is richer 
in H 2 O. The latter could be due to a higher 0:C ratio, sug¬ 
gesting RT Vir had a higher progenitor mass and/or is more 
highly evolved than the Miras (since all 4 stars are so nearby 
that the overall Galactic metallicity gradient has little effect 
although local fluctuations are possible). Cooke & Elitzur 
(1985) showed that, for unsaturated masers, d$/dr increases 
by up to an order of magnitude if the fractional abundance of 
H 2 O, f— [H 2 O] / [H 2 ], increases by a factor of 2. n q is weakly 
dependent on /, so the clouds around RT Vir could be even 
denser by a factor of a few tenths. This could increase the 
acceleration of the gas, but only slightly, as the momentum 
coupling with the dust is already good for most of the den¬ 
sity range in the H 2 O maser clouds (Richards et al. 1999). 

For IK Tau and RT Vir we were able to estimate the 
photon luminosity as a function of distance from the star 
(Section 4.2). This is shown in Fig. 8. Comparing our results 
with the models by Cooke & Elitzur (1985) we found that 
the emissivity was 100-1000 times greater than predicted 
for the mass-loss rates based on CO, Mco . This difference 
can be explained in terms of the clumping of the mass-loss 
into the discrete clouds we observe. The collision rate (and 
hence the maser pump rate) is proportional to p 2 , where p 
is the local density. Thus, the postulated over-density of the 
maser clouds should increase the maser luminosity at the 
same radius. 

From Table 8, the IK Tau maser clouds are over-dense 
by a factor f p ~ 10 — 50, suggesting that the emission at 
an equivalent distance from the star should be > 100 times 
brighter than the value deduced from Mco- This is con¬ 
sistent with the models of Cooke & Elitzur (1985) for a 
fractional H 2 O density of / < 2 x 10 -4 . However, other fac¬ 
tors may be involved as the profile of d<f>/dr as a function 
of r is very irregular, which implies changes in the rate or 
composition of mass loss from IK Tau as well as possible 
axisymmetry of the CSE. 

For RT Vir, the clouds are over-dense by a factor of 
f p ~ 115. This suggests that d^/dr should exceed the value 
deduced from Mco by a factor of roughly 115 2 /(3 2 ) 2 , i.e. 1 
200. In fact, RT Vir is brighter than the models of Cooke & 
Elitzur (1985) by 3 orders of magnitude, even for the highest 
fractional H 2 O density considered, / = 6 x 10~ 4 . 
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Table 8. The estimates of mass loss rate from H2O maser measurements compared with those using CO observations, the maser cloud 
filling factor and other properties. See Section 5.3 for details. 


Source 

(i) 

Vh 2 0 

(10“ 6 M 

(2) 

© yr -1 ) 
(3) 

fp 

(4) 

fv 

(5) 

Mh 2 o 
(10- 9 M 0 ) 
(6) 

tcs 

(yr) 

(7) 

t s 

(yr) 

(8) 

N P 

(9) 

Mtot h 2 0 

(io- 6 m q ) 

(10) 

Mtot inj 

(10- 6 M o ) 

(11) 

IK Tau (a) 

117 

19 

50 

0.0043 

257 

3.0 

26 

11 

66 

68 

IK Tau(b) 

25 

10 

10 

0.0043 

62 

3.0 

29 

10 

16 

75 

U Ori 

25 

12 

110 

0.0014 

141 

2.4 

24 

1 

2 

6 

RT Vir 

15 

2 

115 

0.0044 

44 

1.7 

14 

2 

2 

2 

U Her 

83 

14 

240 

0.0008 

104 

2.0 

23 

1 

4 

8 


6 SUMMARY AND CONCLUSIONS 

MERLIN 22-GHz observations of the 4 low-mass AGB stars 
IK Tau, U Ori, RT Vir and U Her show that the H 2 O masers 
surrounding each star are found in thick expanding shells of 
a few hundred mas in diameter. The inner and outer limits to 
the maser distributions were derived, under the assumption 
of spherically symmetric expansion. The inner radii n are 
at 6-16 au, with outer limits ~4 times greater. The maser 
velocity increases two- or three-fold over this distance, giving 
a logarithmic velocity gradient e in the range 0.5 < e < 1. 

For the first time, the individual 22-GHz masers in each 
CSE have been resolved. The unbeamed H 2 O vapour cloud 
diameter Z is about 2-4 au at a distance of 10-25 au from 
the star. If the clouds are formed close to the stellar surface, 
assuming linear expansion, this corresponds to an original 
cloud radius of ~ 0.1 R,. We find that l is less than the veloc¬ 
ity resonance length for spherical clouds. The average cloud 
velocity span is 0.8-1.7 km s -1 , with a FWHM between 0.5- 
0.8 km s _1 . 

The water vapour cloud density at n is one or two orders 
of magnitude greater than the mean density derived from 
the mass loss rates inferred from CO and other observations. 
The volume filling factor of the shells is < 0.5 per cent. I 11 
addition, we find that the clouds are vulnerable to disruption 
by shocks during the passage through the shell, as the sound 
crossing time (tea) of the clouds and the time the cloud takes 
to travel through the maser shell ( t s ) are related by t cs ~ 
(l/10)t s - These facts suggest that the masers are density 
bounded. We find that the maser survival times and sizes 
of clouds around these AGB stars are ~ 10 times smaller 
than those around RSG, approximately the same ratio as the 
stellar masses and sizes. This implies that it is the properties 
of the central star that determine maser cloud sizes, rather 
than behaviour intrinsic to the wind. 

U Ori and U Her have poorly-filled H 2 O maser shells 
and both have shown OH maser flares. Some individual 22- 
GHz maser components in U Ori have a large beamed size 
which does not decrease with increasing Tb, indicative of 
masing from a slab rather than a sphere. In U Her, the 
masers form a hollow ring. The appearance of the 22-GHz 
masers in both of these sources suggests that they emanate 
from shock-compressed shells, producing strong tangential 
amplification; such shocks could be the cause of the OH 
flares. 

IK Tau and RT Vir have well-filled H 2 O maser shells. 
Observations of both of these sources made ^ 10 years apart 
show a persistent E-W offset between moderately red- and 
blue-shifted emission, although individual masers do not sur¬ 


vive more than ~ 1.5 yr. This can be explained by an equa¬ 
torial density enhancement in a spherical shell, so that the 
brightest masers lie in an oblate spheroid and the equato¬ 
rial plane is at an angle i to the line of sight, such that 
7 t/2 > i > 7r/4. At any one epoch, the shells of U Her and 
U Ori appear asymmetric, but after comparing our maps 
with those published at different epochs (see Section 4 for 
references), we cannot distinguish any systematic deviations 
from poorly-filled spherical shells. 

We are obtaining multi-epoch maps of H 2 O and OH 
masers using MERLIN and the EVN to provide similar res¬ 
olution for all species. This will allow us to study the evolu¬ 
tion of individual maser clouds. We will measure their proper 
motions and distinguish between acceleration and any lati¬ 
tude dependence of expansion velocity e.g. a bipolar outflow. 
We will thus examine the stability of the velocity field in the 
CSE as a whole and maser amplification within individual 
clouds. We will measure changes in maser size and brightness 
with time and so constrain the degree of maser saturation 
and the local overdensity of maser clouds as a function of 
radius. We will develop a fuller model of the velocity fields 
to compare with predictions based on dust properties and 
distribution. 
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